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Abstract 

The fragmentation of b quarks into B mesons is studied with four million hadronic Z 
decays collected by the ALEPH experiment during the years 1991-1995. A semi-exclusive 
reconstruction of -B — iuD^*^ decays is performed, by combining lepton candidates with 
fully reconstructed D^*^ mesons while the neutrino energy is estimated from the missing 
energy of the event. 

The mean value of x^'^, the energy of the weakly-decaying B meson normalised to the 
beam energy, is found to be 

{x'^'^) = 0.716 ± 0.006 (stat) ± 0.006 (syst) 

using a model-independent method; the corresponding value for the energy of the leading 
B meson is {x^) = 0.736 ± 0.006 (stat) ± 0.006 (syst). The reconstructed spectra are 
compared with different fragmentation models. 
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1 Introduction 



The process of hadron production at e^e~ colliders is usually modelled as the convolution 
of a perturbative part (hard gluon radiation for energies above approximately 1 GeV) 
and a non-perturbative part, called hadronisation or fragmentation, in which the quarks 
are confined in colourless hadrons. While the first step is in principle calculable, the 
fragmentation needs a phenomenological approach and is usually parametrised in terms 
of the variable 



hadron 



1 



quark 



where p\\ is the hadron's momentum along the direction of the quark, and {E + p)quark 
the sum of the quark energy and momentum just before fragmentation, i.e. taking into 
account initial and final state photon radiation, and hard gluon emission. 

With this definition, the fragmentation process can be described in terms of the prob- 
ability of a hadron H to be generated with a given z, called D^[z), where q is the flavour 
of the generating quark. In this paper the fragmentation of b quarks is studied. 

The fraction z is not accessible experimentally, and hence a direct reconstruction of 
D^[z) is not possible. The energy spectrum of b hadrons is therefore described in terms of 
the scaled energy x^, defined as the ratio of the heavy hadron energy to the beam energy 



Xb = 



had 



(2) 



beam 



In contrast to the z variable, the effects of initial and final state radiation and hard gluon 
emission are not unfolded. 

In the analysis presented, the energy of B mesons is reconstructed using a partially 
exclusive method: semileptonic decays B iuD^*^ are identified by pairing lepton can- 
didates with fully reconstructed D^*^ mesons; the scaled energy of the weakly-decaying 
B meson is then computed adding an estimate of the neutrino energy. Five channels are 
chosen because of their good signal purity and statistical significance; they are shown in 
Table 



In the following x°^'^ indicates the reconstructed energy of B meson candidates. 



X 



wd 
B 



the energy of weakly decaying B mesons, corrected for detector acceptance and resolution; 
x^ stands for the corrected scaled energy of the leading B meson, that is the first meson 
produced in the fragmentation process, which can also be a heavier resonance {B*, B**). 



Channel 
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B'- 
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D'- 


Kn 
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B°- 


iuD* 


D* - 




D^- 


Kmrn 
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B - 


> iuD^ 






D^- 


Kn 
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B°- 


£uD 






D - 


> Kmr 
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B^- 


iuD* 


D" - 




D^- 





Table 1: i?-decay channels used in the analysis. 



The analysis uses the full LEP I statistics collected by ALEPH between 1991 and 
1995, amounting to almost four million hadronic Z decays. Recently this data set has 
been reprocessed using improved reconstruction algorithms. The main benefits for this 



analysis are related to the enhanced secondary vertex reconstruction efficiency and the 
improved particle identification. A discussion of the reprocessing can be found in |Q. 

After a description of the ALEPH detector, the selection of i? — > £uD^*^ decays is 
detailed in Section ^. In Section ^ the reconstruction of the B meson energy is described, 
followed by the extraction of the spectrum and comparison with the predictions of different 
models in Section |[ Systematic errors are discussed in Section |^, and checks on the self- 
consistency and robustness of the analysis are presented in Section ^ 

2 The ALEPH detector 

The ALEPH detector and its performance are described in detail elsewhere p|, ||. A high 
resolution vertex detector (VDET) consisting of two layers of silicon with double-sided 
readout measures r0 and z coordinates at average radii of 6.5 cm and 11.3 cm, with 
12 yum resolution at normal incidence. The VDET provides full azimuthal coverage, and 
polar angle coverage to |cos6'| < 0.85 for the inner layer and |cos6'| < 0.69 for both 
layers. Outside VDET, particles traverse the inner tracking chamber (ITC) and the time 
projection chamber (TPC). The ITC is a cylindrical drift chamber with eight axial wire 
layers with radii between 16 and 26 cm. The TPC measures up to 21 space points per track 
at radii between 40 and 171 cm, and also provides a measurement of the specific ionization 
energy loss (dE/dx) of each charged track. These three detectors form the tracking 
system, which is immersed in a 1.5 T axial magnetic field provided by a super-conducting 
solenoid. The combined tracking system yields a transverse-momentum resolution of 
<j{Pt)/pt = 6 X 10~^pT{GeV / c) © 0.005. The resolution of the three-dimensional impact 
parameter for tracks having two VDET hits can be parametrised as a = 25 /im + 95/im/p, 
{p in GeV/c). 

The electromagnetic calorimeter (ECAL) is a lead/wire chamber sandwich operated 
in proportional mode. The calorimeter is read out in projective towers that subtend 
typically 0.9° x 0.9°, segmented in three longitudinal sections. The hadron calorimeter 
(HCAL) uses the iron return yoke as absorber. Hadronic showers are sampled by 23 
planes of streamer tubes, with analogue projective tower and digital hit pattern readouts. 
The HCAL is used in combination with two double layers of muon chambers outside the 
magnet for muon identification. 

3 Selection of S ^ ii^D^*^ decays 

A Monte Carlo simulation based on JETSET 7.4 Q and tuned to ALEPH data Q has 
been used in order to extract resolution functions, acceptance corrections and background 
compositions. About five million bb events were simulated, and more than twice the data 
statistics of qq events. The present analysis uses bb events to determine the x^'^ and 

spectrum starting from observed spectra, and qq events to evaluate the non-66 
component of the selected sample. 

The decays B iuD^*^ are searched for in hadronic events, containing at least one 
lepton (electron or muon) identified using standard criteria 0. The momentum cut used 
to define lepton candidates is p > 2 GeV/c for electrons and p > 2.5 GeV/c for muons. 
The transverse momentum pj- of the lepton with respect to the nearest jet, with the lepton 
excluded from the jet, is required to be larger than 1 GeV/c, which helps rejecting fake 
candidates and leptons not coming from direct decays of b hadrons. Both electron and 
muon candidates are required to have a measured dE/dx compatible with the expected 



value. 

Events are divided into two hemispheres by a plane perpendicular to the thrust axis; 
in each hemisphere containing a lepton a D meson reconstruction is attempted in the 
decay modes described in Table |I|. At least two charged tracks from the D meson decay 
are required to have VDET hits, in order to ensure a good reconstruction of the D 
vertex position and to reject combinatorial background. Loose cuts are applied to track 
momenta, in order to minimise the bias in the B momentum distribution. Tracks are not 
considered as kaon candidates if their measured ionization is incompatible with the kaon 
hypothesis by more than three standard deviations. The charge of the kaon candidate 
is required to be the same as that of the lepton, as expected for semileptonic B meson 
decays. 

Tracks assigned to a D meson decay are fitted to a common vertex, and the track 
combination is rejected if the of the fit is larger than 20. If more than one combination 
fulfils this requirement for channels 3 and 4, the one with the smallest is chosen. In 
channel 5, the vr" closest in angle to the charged pion is selected and added to form the 

For channels 1, 2 and 5, a soft pion tTs is added to the D candidate to form a D* meson; 
the TTs momentum is required to be larger than 250 MeV/c and smaller than 3 GeV/c. The 
difference between the reconstructed D* and D masses is required to be within 5 MeV/c^ 
of its nominal value. In the case of multiple candidates in a given hemisphere, the track 
combination is chosen for which the reconstructed {D* — D) mass difference is closest to 
the nominal value. 

A vertex fit is performed using the D candidate and the lepton track, and again the 
combination of tracks is rejected if the of the fit is larger than 20. The B vertex is 
required to lie between the interaction point, reconstructed event-by-event, and the D 
vertex. 

Channels 3 and 4 are further enriched in signal events using harder cuts on the kaon; in 
addition a tTs veto is applied: if a track is found which is compatible with the reconstructed 
B vertex and the combination track-D candidate has a mass close to the mass of the D*, 
the candidate is discarded. This procedure reduces the overlap between channels at the 
permil level. Finally, tighter cuts on the reconstructed D mass and the of the vertex 
fit are imposed. 

The D mass spectra are shown in Fig. ^ The reconstructed D mass peaks are fitted 
in a region between 1.7 and 2.0 GeV/c^ with a Gaussian and a linear component. Table |] 
shows the chosen D mass windows, the number of reconstructed candidates and the fitted 
Gaussian fractions. 

The fractions of the Gaussian components measured in the Monte Carlo are compatible 



Channel 


D window 
(MeV/c2) 


Events 


Resolution 

(MeV/c^) 


Gaussian 
Fraction 


D* - 


DOyr, - 


^ Kn 


1864 ±30 
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8.3 


89 % 


D* - 


^Ovr, DO - 




1864 ±30 


388 


6.2 


69 % 






^ Ktt 


1864 ± 15 


1079 


8.4 


81 % 




D 


Ktitt 


1869 ± 30 


580 


7.4 


64 % 


- 


- 




1864 ±50 


693 


25 


63 % 



Table 2: For the five channels, the D mass window, the number of events in the window, the mass 
resolution and the Gaussian fraction. 
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Figure 1: Reconstructed D mass peaks in the five channels. In the second channel the second peak at 
lower mass comes from the decay channel Z?" Kinnnr'^ and is excluded from the fit. 



with those in the data within statistical errors, while the widths are about 5 — 10% smaller; 
this is taken into account in the evaluation of the systematic uncertainties. 

4 B energy reconstruction 

The scaled energy of the weakly-decaying B — > iuD^*^ hadron is estimated as 

-^beam 

The terms and are provided by the direct reconstruction, while the neutrino 

energy E^ is estimated from the missing energy in the hemisphere: 

TP — fi^hemi _ 7-.hemi fA\ 

~ -^tot -^VIS 5 \^) 

where -E^™' is estimated taking into account the measured mass in both hemispheres |^: 

2 2 
TTihemi _ 771 1 same ""oppo 

-^tot - -^beam H TT^ ' ' ' 

'^-^-'beam 



Both charged and neutral particles are used in Eqns. (Q) and (^. In the lepton hemisphere 
neutral hadronic energy is expected to come only from fragmentation. Therefore, in order 
to avoid spurious calorimetric fluctuations, its contribution is taken into account only 
outside a cone of 10 degrees of half opening angle around each of the B meson decay 
products. Table ^ shows the resolution on xb estimated on simulated bb events; the 
distributions are well described by two Gaussians, accounting for core and tails. 



Channel 


Core (%) 
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resolution 


Tail 
resolution 
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Table 3: For the five channels the xb resolution on simulated events. The resolution can be parametrised 
with two Gaussians, describing the core and the tails. 



5 Unfolding methods 

The scaled energy of the weakly-decaying B mesons, x^^, is reconstructed in 20 bins 
between and 1 with a variable width. In each bin, the non-66 background is estimated 
using the simulation, and subtracted from the spectrum. This amounts to about 2% of 
the events, concentrated mostly at low x^*^. The measured spectra after subtraction are 
shown in Fig. 0. 

With these events two different kinds of analyses can be performed: 

• a model-dependent analysis, in which different fragmentation models available in 
the literature are tuned to fit the observed spectra; 

• a model-independent analysis, in which the shapes of x^'^ and x^ are reconstructed 
by correcting the observed spectra for detector acceptance, resolution and missing 
particles. 



5.1 Model-dependent analysis 

Various fragmentation functions D^{z) are implemented in the Monte Carlo generator 
JETSET 7.4, which also simulates initial and final state photon radiation and hard gluon 
emission. The reconstructed spectra obtained from the simulated 66 events are tuned to 
best reproduce the x^^ distribution observed in the data, by minimising the global x^- 
The following parametrisations for D^{z) are used: 

'2 



Peterson et al. [9] : 
Kartvelishvili et al. ||lO| : 
Collins et al. [ITT 
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Figure 2: In the five channels, x'^^ spectra, after non-66 background subtraction and before acceptance 
corrections. 



The minimisation is performed with respect to the free parameter of each model, and 
the is written as: 



X 



5 20 

EE 

C=l 4 = 1 



(6) 



where c runs over the channels used, i runs over the x^'^ bins defined as in Table | 



n 



DT 



and n'^*^ are the number of candidates per channel and per bin observed in the data and 
expected from the Monte Carlo, normalised to the same number of entries. The quantities 
a are defined as statistical uncertainties. 

Table ^ shows the fitted values for the different model parameters, together with 
statistical and systematic uncertainties from the sources discussed in Section ||. Also 
shown are the values for the mean scaled energy. 



Model 


Fit results 


Mean energies 


Peterson 
i 


= (31±3±5)xl0-^ 
XVA^DOF = 117/94 


= (700 ± 4 ± 5) X 10"^ 
(xL) = (721±4±5)xl0-3 


Kartvelishvili 
10 


dp = 13.7 ±0.7 ±1.1 
XV^^DOF = 107/94 


= (713±4±6)xl0-^ 
(x^) = (734±4±6)xl0-3 


Collins 

|11| 


eb= (185 ± 25 ± 41) X 10-^ 
xVA^dof = 181/94 


{xl"") = (681 ± 4 ± 5) X 10"^ 
(x^) = (701 ±4±5)xl0-3 



Table 4: Fit results with different fragmentation models. The systematic errors account for the sources 
of uncertainties discussed in Section ^. The x^/Noof is calculated using statistical errors only. 



The Kartvelishvili model describes the data slightly better than the Peterson model. 
The Collins model is clearly disfavoured. 



5.2 Model-independent analysis 

The xl'^ and spectra are obtained by correcting the observed x'^'^ spectra for accep- 
tance, detector resolution and missing particles. 

The normalised binned spectrum /j (^x^^^^ can be obtained using the relation 

1 5 1 20 

/^(-S^) = ^E-^EG-(c)nnc) , (7) 

where ef^i^c) is the acceptance correction in bin i for channel c; n^"'"(c) is the number of 
reconstructed B mesons in the data for channel c, with a measured energy falling in bin 
j; GJ^^{c) is the resolution matrix that links mesons with x'^^ in bin j and x^'^ in bin i, 
for channel c; T is the normalisation factor defined by the condition J2i fi = ^- A similar 
equation holds for the extraction of (x^^ , where the effect of the missing particles from 
B* and B** decays is folded in and G^. 

The acceptance corrections and the resolution matrix Gij are taken from the simula- 
tion. The acceptance corrections show a different behaviour among the different channels; 
the two extreme situations are shown in Fig. ^. A dependence on the fragmentation func- 
tion present in Monte Carlo is induced in the measured spectrum through Gij] hence, the 
Monte Carlo used to calculate G^ must be reweighted to the best estimate of fi (x^*^^ 

from data. This is done using an iterative procedure, calculating //^(x^'^) using the G^"^ 
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Fi gure 3: Acceptance corrections ej^'^(c) for x^*^. The absolute scale is chosen as to conserve the total 
number of selected events for each channel. Only the distributions for channels 1 and 4 are shown, since 
they represent the extreme behaviours. 



from Monte Carlo reweie hted to /.^-^[x'S^). The weights Wi = //^ ^(DT)//j(MC) are ap- 
phed to standard Monte Carlo events; to avoid fluctuations due to the limited statistics 
in data events, the distribution of the weights Wi is smoothed with a polynomial function. 
Possible systematic effects related to the smoothing are studied in Section ^ 

The whole procedure is then repeated until the change in /j in consecutive iterations 
is a small fraction of the statistical errors. 

The statistical error matrix Eij is calculated by repeating 20 x 5 analyses, varying in 
each of them the quantities ?t,^"'"(c) by one standard deviation: 

E^,= E E (/;^^'^-/f™)(/r-/f°) , (8) 

c=l,5 fc=l,20 

where fj'^^'^ is the result of the convergence for fi when n^^(c) is varied by its statistical 
error, and ff™ the nominal result. 

The results, together with the statistical and systematic errors, are given in Table |^. 
The full error matrices are shown in the Appendix. 

From the binned spectra fi{x'^'^^ the mean value is calculated as 

20 

(xb') = E^^/«(^5') , (9) 

i=l 

where Xi is the central value of bin i. The bin size chosen is such that the deviation from 
linearity of the distribution within a bin is negligible. 

The statistical error on (xb) is calculated using the same procedure used as for fi. 

The results for (x^'^) and (x^) are 

(xS'^) = 0.7163 ± 0.0061 (stat) , 
(x^) = 0.7361 ± 0.0063 (stat) . 



Process 


BR(%) 


B Diu 


1.95 ±0.27 


B DHv 


5.05 ±0.25 


B DWX£z/ 


2.7 ±0.7 


with B Diiu 


0.63 ±0.11 


with B D^iu 


0.23 ±0.09 


b —>■ uiu 


0.15 ±0.10 




9.85 ±0.80 


Inclusive B — > ipX 


10.18 ±0.39 



Table 5: Exclusive branching ratios for the B IvX process ||T2|, The sum is consistent with the 
measurement of the inclusive B ivX rate. 



6 Systematic errors 

Possible systematic effects due to uncertainties on the physics parameters used in the 
Monte Carlo, limited accuracy in the simulation of the detector response, or effects in- 
trinsic to the analysis method have been investigated. 

The physics parameters used in the Monte Carlo simulation that are relevant for the 
analysis are adjusted to the most recent experimental measurements and varied within 
their estimated uncertainty by reweighting simulated events. The effect of the reweight- 
ing propagates to the results through the resolution matrix Gjj(c) and the acceptance 
corrections ej(c), which are taken from the simulation. The differences from the standard 
results are taken as systematic errors. 

The sources of uncertainty considered are: 

• Semileptonic decays of B mesons. 



The current experimental knowledge [jT2|, [T^ of the semileptonic branching ratios of 



B mesons is summarised in Table ^. The sum of the exclusive (or semi-exclusive) 
rates is consistent within errors with the inclusive measurement of BR(i? — > IvX). 

The analysis is not sensitive to the total BR(i? — > IvX), but is affected by a change 
in the relative rates of the different components, since these contribute in a different 
way to the average acceptance corrections and resolution matrix. 

Six sources of systematic error are calculated using the values in Table ^: 
1. The inclusive BR(i? D^*^Xtv) is varied within its experimental error. 

A(a;;^'^) = 0.0019 A(x^) = 0.0020 



2. The rate for the narrow Di state is varied by its experimental error, while 
leaving the total BR(i? D^*^X£u) at its central value. 

Aix'S^) = 0.0001 A(x^) = 0.0001 

3. The rate for the narrow D2 state is varied by its experimental error, while 
leaving the total BR(i? D^^^Xiu) at its central value. 

Aix'f) = 0.0001 A(a;^) = 0.0001 



4. The rate of wide D** states, not yet measured, is put to zero and compensated 
with non-resonant D^*Hutc final states, thus leaving the total BR(i? — > D^*^Xiu) 
at its central value. 

A(x;^'^) = 0.0017 A(x^) = 0.0016 



5. The BK{B Diu) is varied by its experimental error: 

Aix'f) = 0.0008 A(x^) = 0.0007 

6. The BR(i? D^iu) is varied by its experimental error: 

A{x'^'^) = 0.0008 A{x^) = 0.0008 



Missing particles from B** production. 

When deriving the energy spectrum of the leading B meson, the correction due the 
energy carried away by the pion produced in the B** decay enters in the resolution 
matrix and the acceptance corrections. The rate of 6 — B** is varied within its 
experimental error: fs** = 0.299 ± 0.058 [H, 15, 1^, and the resulting systematic 



error is A(xg) = 0.0025. The weakly-decaying B meson spectrum is not affected 
by this source of uncertainty. 

Modelling of the B** production. 

From spin counting, the relative rates of (-Bi, Bq, B*, B2) are predicted to be 
(3,1,3,5) [1^; changing this to (1,1,1,1) gives a systematic error of A(a;^) = 0.0004. 

Production of B* from b quarks. 

Due to the small mass difference between B* and B mesons, the effect of B* pro- 
duction in b quark fragmentation is found to be much smaller than for B**, and it 
is completely negligible for the present analysis. 



The relevant sources of systematic uncertainties due to the detector simulation are 
identified to be: 

• Neutrino energy reconstruction. 

The accuracy of the neutrino energy reconstruction is checked in hadronic events 
enriched in light primary quarks and in hadronic decays of rr events. In the first 
sample, a "fake" neutrino is simulated by removing a charged particle from the 
reconstructed event; its energy is then reconstructed using Eqns. (|) and (^. The 
method can be applied both to data and Monte Carlo events, determining the bias 
between the reconstructed energy and the momentum measured with the tracking 
system. Such a bias is found to be reproduced by the Monte Carlo with a precision 
better than 50 MeV, for all momenta of the deleted track. In the second sample, 
where the event topology is much simpler, the energy of the "reconstructed" Ur 
is compared between data and Monte Carlo events. Also in this case the worst 
discrepancy observed is smaller than 50 MeV. This value is used as a conservative 
estimate of the systematic uncertainty on the neutrino energy, resulting in A(a;^*^) = 
0.0023 and A(x^) = 0.0023. 



Vertexing and charm meson reconstruction. 

If the purity and the kinematic properties of the selected candidates are not well 
described by the simulation, the acceptance corrections and resolution matrices 
can be inadequate. In order to check for these effects, the distributions of the 
probability for the reconstructed D vertices are compared, channel by channel, with 
the simulation. Small differences are observed, and the Monte Carlo distribution is 
reweighted in order to reproduce the data. The shift in the corrected average energy 
is taken as systematic uncertainty. The resulting error estimates are A(x^'^) = 
0.0001 and A(x^) = 0.0001. 

Furthermore, the reconstructed D mass distributions in data and Monte Carlo are 
compared. In simulated events the widths of the mass spectra are found to be 
5 — 10% smaller, while the fractions of the Gaussian components, estimated from a 
fit to the sidebands, are reproduced within their statistical error of about 5%. The 
mass cuts reported in Table |]are adjusted in order to take into account both effects, 
taking the total shift in the extracted energy spectrum as systematic uncertainty. 
The resulting estimates are 



width = 0-0011 



/\{xl^) = 0.0021 A(x^) 



L ' 



, = 0.0010 

width 

. = 0.0021 

purity 



purity 

Possible systematic effects related to the analysis procedure are: 

• Background subtraction. 

As previously explained, a bin-by-bin subtraction of candidates not coming from 
bb events is performed before deriving the B meson energy spectra. The efficiency 
for this kind of background has been extracted directly from data events, which 
have been enriched in background events by selecting wrong sign candidates. It is 
found to be compatible with Monte Carlo simulation within the statistical error of 
about 25%. The background is varied within this range and the systematic errors 
associated are A(a;;g'^) = 0.0021 and A(x^) = 0.0022. 

• Monte Carlo statistics. 

Statistics of simulated events are larger than for data events by a factor of 5. In 
order to evaluate the related uncertainty, the acceptance corrections ei(c) and the 
matrix elements Gjj(c) are varied randomly by their statistical error in a series of toy 
experiments. The scatter of the results for (x^'^) and (x^) is taken as an estimate of 
the uncertainty due to the limited Monte Carlo statistics, yielding A(x^'^) = 0.0029 
and A(x^) = 0.0031. 

Adding in quadrature all the systematic contributions, the final results are: 

{x'^'^) = 0.7163 ± 0.0061 (stat) ± 0.0056 (syst) , 

(x^) = 0.7361 ± 0.0063 (stat) ± 0.0063 (syst) . 

The bin-by-bin results for the measured spectra, with the total systematic uncertain- 
ties, are shown in Table |^, while the statistical and total error matrices are reported in 
the Appendix. The spectra are also shown in Fig. H, where they are compared with the 
Monte Carlo predictions from different fragmentation models, with the free parameters 
fitted to the data. 

The models of Peterson and Kartvelishvili give the best agreement with the data, and 
are compared with the x^*^ measurement in Fig. |^. 
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Figure 4: Scaled energy of the leading and weakly-decaying B meson, as reconstructed from data. 

The best-fit distributions for the Peterson model, the Kartvelishvili model, and the Collins model are 
superimposed. For the data, the bin-to-bin errors are highly correlated, as shown in the error matrices 
in the Appendix. 
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Figure 5: Scaled energy of the weakly-decaying B meson, as reconstriictcd from data. The inner error 
bars represent statistical errors, the larger ones the total uncertainties. The best-fit distributions for the 
Peterson model and the Kartvelishvili model are superimposed. For the data, the bin-to-bin errors are 
highly correlated, as shown in the error matrices in the Appendix. 
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Table 6: The extracted spectra for the weakly-decaying B meson and the leading B meson. All numbers 
are given in units of 10^"^. 

7 Systematic checks 

Possible systematic effects intrinsic to the analysis method are checked by measuring the 
energy spectra in a sample of 8 million simulated qq events. The average values for the 
scaled energies of the weakly-decaying and leading B mesons are measured to be: 



(x^)mc = 0.711 ± 0.005 (stat) , 
{x'1^)mc = 0.692 ±0.005 (stat) , 

which compare well with the true values: 

(4)m'S = 0.712 , 
{^t)MC = 0.692 . 

Electron and muon identification are affected by different sources of background, and 
the selection efficiencies and purities have a different dependence upon the track kine- 
matics and isolation. It is therefore interesting to perform the analysis using separately 
events with electron candidates or muon candidates. Consistent results within uncorre- 
lated errors are found: 

W)electrons = 0.724 ± 0.010 , (x^ ) electron. = 0.743 ± 0.010 , , . 

(x5<^)„s = 0.700 ± 0.014 , (x^)„ = 0.720 ± 0.014 . ^ > 

The acceptance corrections for the five channels are significantly different (Fig. ^ and 
the same is true for the resolution matrices. An inaccurate description of these inputs 



would easily lead to incompatible results among the different channels. This is checked 
by performing the analysis separately in the five sub-samples. The results, reported in 
Table are compatible within uncorrelated uncertainties. 



Channel 




(^^^) 


1 


0.700 ±0.015 


0.720 ±0.016 


2 


0.700 ±0.020 


0.720 ±0.022 


3 


0.714 ±0.012 


0.733 ±0.013 


4 


0.720 ±0.019 


0.740 ±0.019 


5 


0.738 ±0.012 


0.755 ±0.013 



Table 7: Results using the five channels separately. The errors are uncorrelated. 



It has been checked that the results are independent of the choice of fragmentation 
functions in the Monte Carlo sample used to estimate the resolution matrix Gij and the 
acceptances e^. 



As explained in Section |5]^, the weights applied to reweight Gij to a given frag- 
mentation function are smoothed with a polynomial function to reduce the bin-to-bin 
fluctuations. However, the values for the mean scaled energies move by a small fraction of 
the statistical errors when such smoothing is not applied, and the total statistical errors 
remain nearly constant: 



)= 0.7177 ±0.0060 , (x^) = 0.7370 ±0.0065 . (11) 



Heavy flavoured hadrons originating from gluon splitting g ^ bb have an energy much 
lower than hadrons coming from primary b quarks. A check on Monte Carlo events shows 
that the contribution of such events is negligible. 

The analysis uses a binned representation of the fragmentation functions to compen- 
sate the relatively small statistical sample in the data. The binning chosen must not 
introduce biases in the measured values nor should it affect the statistical errors. This 
is checked by performing a number of analyses in which the binning is varied randomly 
around the standard one. Both the central values and the statistical uncertainties are 
stable. 



8 Conclusions 

Using the data collected by the ALEPH experiment at and around the Z resonance in 
the years 1991-1995, about 3400 semileptonic and decays have been selected. The 
scaled energy spectra of weakly-decaying and leading B mesons have been reconstructed, 
and their mean values were found to be: 



(4) 



0.7163 ± 0.0061 (stat) ± 0.0056 (syst) 
0.7361 ± 0.0063 (stat) ± 0.0063 (syst) 



The observed spectra have been compared with the prediction of JETSET 7.4 using 
different fragmentation models. The models of Peterson et al. |^ and Kartvelishvili et 
al. |]TU[ give a reasonable description of the data while the Collins et al. 
clearly disfavoured. 



nfl model is 



This measurement supersedes a previous analysis from ALEPH [IS], which used a 
different method and smaller statistics. 

The present result is compatible with the published results using b hadrons from 
L3 ig, OPAL [gg and SLD 0], and using B mesons from OPAL [0. 
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A Appendix 

Tables |^ and ^ show the statistical and total error matrices for x^'^; Tables ^ and [Tl] give 
the same information for x^. 
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Table 8: Statistical error matrix for x'^'^. All the numbers are in units of 10 ^. 
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Table 9: Total error matrix for x'g^. All the numbers are in units of 10 ^. 
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Table 10: Statistical error matrix for x^. All the numbers are in units of 10 
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Table 11: Total error matrix for x\. All the numbers are in units of 10 ^. 
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